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Title 

Frequency Estimation 

Technical Field 

5 This invention concerns a method for estimating the frequency of a 

single frequency complex exponential tone in additive Gaussian noise. In 
another aspect the invention is a frequency estimation program for estimating 
the frequency of a single frequency complex exponential tone in additive 
Gaussian noise. In further aspects, the invention is computer hardware 
10 programmed to perform the method. 

Background Art 

Earlier work to estimate the frequency of a single frequency complex 
exponential tone in additive Gaussian noise uses the fast Fourier transform 
15 (FFT) algorithm. The initial work on this topic was introduced by Rife and 
Boorstyn [1-3]. This paper introduces an algorithm employing the FFT, which 
produces an estimate of the frequency with extremely low variance of the error. 
The variance of the frequency estimate is independent of the frequency of the 
signal. The algorithm has a low computational complexity implementation. 

20 

The received signal, r[n], is given by, 

r[n] = s[n] + T)[n], for n = 0,1,2 N-1 (1) 

where: 

s[n] = Ae J ** T ' 1 

i 7 ? M}** 1 is a set of independent, complex, zero mean, Gaussian random 

25 

variables with variance a 2 , 
7 R [n]=Re{y 7 [n]} > 

77 1 [n] = Imag{7[n]}, 

f is the frequency of the tone, 



30 T s is the sampling period, 



-y = var[7 7R [n]]=var[7 I [n]] 



and, A is the signal amplitude. 



The sampling frequency, f s , is given by, 
5 f s samples/s. (2) 

The signal to noise ratio of each complex signal plus noise sample is given by, 



SNR=— (3) 



A 2 
a 

Rife and Boorstyn [1-4] suggest a method of estimating f by using a FFT. 
10 It is assumed that 0 < f < f s . First, a coarse search is performed. Under 
noiseless conditions, the absolute value of the FFT output coefficient 
corresponding to the bin centre frequency closest to f will be maximum over the 
set of absolute values of the FFT output coefficients. The coarse search, 
performed by the FFT, narrows the frequency uncertainty, to 

15 Hz, where an N point FFT is used. Then, a fine search method is used to 

further reduce the frequency uncertainty. A secant method is used to compute 
the estimate of f by successful approximates. 



Define, 





>(0) ' 




~Y(0) " 




r(1) 




Y(1) 


20 

r = 


r(2) 


, Y = 


Y(2) 




r(N-1) 




. Y(N-1) 



where, Y =FFT(r) and FFT(. ) is the Fast Fourier Transform Operator. 



25 Then the Rife and Boorstyn coarse search is, 



k^max' {|Y(k)|:0<k<N-l} 

0 Sk SN> 1 



(5) 



and, 

5 fo=(^ L )f« (6) 
where, 

f 0 is the coarse frequency estimate in Hz. 

10 

Numerous other frequency estimation approaches have been suggested 
in the literature [5 -1 0]. 

Disclosure of the Invention 

15 In a first aspect the invention is a method for estimating the frequency of 

a single frequency complex exponential tone in additive Gaussian noise, 
comprising the steps of: 

performing the fast Fourier transform (FFT) on the tone; 
estimating the frequency as the frequency corresponding to the largest 
20 FFT output coefficient magnitude; 

computing a discriminant which is proportional to the frequency error in 
the initial frequency estimate using modified coefficients of the discrete Fourier 
transform (DFT) with center frequencies plus one half and minus one half of the 
FFT bin spacing relative to the original frequency estimate; 
25 mapping the value of the discriminant into the estimate of the frequency 

error in the initial frequency estimate using a mathematically derived function; 

adding the estimated frequency error estimate to the initial FFT estimate 
to get a first interpolated frequency estimate; 

computing a further discriminant which is proportional to the frequency 
30 error in the first interpolated frequency estimate using modified coefficients of 
the discrete Fourier transform (DFT) with center frequencies plus one half and 
minus one half of the FFT bin spacing relative to the first interpolated frequency 
estimate; 





5 



mapping the value of the further discriminant into the estimate of the 
frequency error in the first interpolated frequency estimate using the 
mathematically derived function; and 

adding the estimated frequency error in the first interpolated frequency 
5 estimate to the first interpolated frequency estimate to get a second 
interpolated frequency estimate. 

The first interpolated frequency estimate is quite accurate because it is in 
a region of relatively low noise induced frequency error. The method generates 
an unbiased, low error variance estimate of the frequency. The performance of 
10 the method, above the signal to noise ratio threshold, is about 0.06 dB above 
the Cramer-Rao lower bound. The method is ideally suited to be utilised in a 
number of communications, signal processing and biomedical applications. 
The method is easily implemented in hardware or software with low 
computational overhead. 
15 In theory, this technique of iteratively deriving an interpolated frequency 

estimate and then, using the frequency discriminant, a more precise frequency 
estimate can be continued infinitely times until a fixed point (or solution occurs). 
At this fixed point, the discriminant function has .zero value. 



Several functions have been identified to compute the discriminant. In 



20 practice, different functions may require a different number of iterations to 
essentially converge to a fixed-point solution. However, discriminant functions 
defined by a wide class of functions using two DFT coefficients as the input 
converge to the same solution and therefore exhibit identical noise 
performance. 

25 A first example of the discriminant, or distance metric, of frequency 

estimation error is: 




(10) 



(9) 



30 



and, 



max 



A k 

So for the initial frequency estimate using the FFT, f 0 T s = — —and s = 0. 

N 

Other examples of the discriminant having the properties required for the 
algorithm include: 

B D= Hir-i«r >fory>0-i 
r\fir+\<*r 

and in particular, 

l |/?r--|gp 
2\/}\ 2 +\ar 
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and 



P + a 



15 where Re[.] is the real part and * denotes the complex conjugate. 

In these equations, p and a are the modified DFT coefficients defined by, 

/? = Y(k mw+ |) = 2;r(n)e j ^-^ 



20 



n=0 



1 £J -j2*n- 



and f a = Y(k_-^) = £r(n)e' » 

It is also possible to define discriminant functions which use more than two DFT 
coefficients to obtain further improvements in frequency estimation 
performance in additive Gaussian noise relative to discriminants that use only 

25 two DFT coefficients. An example where 2M+2 coefficients are used, where 

TV 

0zM< — ~l and the FFT coefficients are used in the discriminant with 
2 

optimal weighting coefficients obtained by using the concept of matched 
filtering is, . 



ZC , {Y[(k^+I + m)modlfl-r[(k^-i-m)niodN]} 

j,^^^ 1 ] 

£{C , {Y[(k M +i + m)niodN]+Y # [(k nttX ---ni)modN]} 

m=0 -■*-<"] 2 2 



N 

where, 0 £M £ 1 , mod N indicates modulo N, 

2 



5 and, where, " denotes complex conjugate. 



C = e sm[4 + m)] 

sin[-(i + m)] 

1 0 and, Y^,,^ + i + m) and Y(k nlax - i - m) are the modified DFT coefficients given 
by, 



15 end, Y(k 01K -~m) = 2r(n)e" 

* n=0 



The discriminant using more than two DFT coefficients may be used in the last 
20 iteration to obtain additional frequency accuracy. In a similar manner, 
discriminant functions may be formulated which use more than two DFT 
coefficients and less or equal to all N FFT coefficients. 

In a second aspect, the invention is a frequency estimation program for 
25 estimating the frequency of a single frequency complex exponential tone in 
additive Gaussian noise, wherein the frequency estimation program has 
functionality to perform the method. 

In a third aspect, the invention is computer hardware programmed to 
perform the method. The hardware may comprise a DSP processor chip, or 
30 any other programmed hardware. 
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Brief Description of the Drawings 

Examples of the invention will now be described with reference to the 
accompanying drawings, in which: 

Fig. 1 is a graph that illustrates the FFT Coefficients, where the signal 
5 frequency is closer to the lower FFT frequency than the higher FFT frequency; 

Fig. 2 is a graph that illustrates the FFT Coefficients, there are two equal 
peak coefficients and the signal frequency is halfway between; 

Fig. 3 is a graph that illustrates the FFT Coefficients, where the signal 
frequency is closer to the upper FFT frequency than the lower FFT frequency; 
10 Fig. 4 is a Flow Diagram for the Frequency Determination Algorithm; 

Fig. 5 is a graph that illustrates the ratio of the variance of the normalized 
frequency error, s-s , to Cramer-Rao Bound variance in dB as a function of 
the FFT length, N; and 

Fig. 6 is a graph that illustrates the variance of the normalised estimator 
15 frequency error estimate against the frequency error for the first interpolation. 

Simulations of the invention show the rms frequency error performance 
of the algorithm vs SNR in dB, for different values of N. 

Figures 7-12 include curves for one interpolation, two interpolations, and 
the Cramer-Rao Bound, where: 
20 Fig. 7 is a graph showing RMS normalised frequency error vs SNR for 

N=2; 

Fig. 8 is a graph showing RMS normalised frequency error vs SNR for 

N=4; 

Fig. 9 is a graph showing RMS normalised frequency error vs SNR for 

25 N=16; 

Fig. 10 is a graph showing RMS normalised frequency error vs SNR for 

N=64; 

Fig. 11 is a graph showing RMS normalised frequency error vs SNR for 
N=256; and 

30 Fig. 12 is a graph showing RMS normalised frequency error vs SNR for 

N=1024. 

* The two-interpolation case essentially achieves the performance of the 
infinite interpolation case. 



35 Best Modes of the Invention 

Referring first to Figs. 1 to 4, the received signal r[n] is given by: 
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r[n) = s[n] + nW, for n = 0,1,2, .... N-1 (1) 

s[n] = Ae i2 ^ T ', 
where: 

{n W)T is a set of independent, complex, zero mean, Gaussian random 
5 variables, 

/7 R [n] = Re{77[n]}, 

^[n^ImagtoM}, 

f is the frequency of the tone, 
T s is the sampling period, 

-y- = var[7 R [n]] = var[7 I [n]] 

10 

and, A is the signal amplitude. 

A fast Fourier transform is performed and the sampling frequency, f 8l is given 
by, 

15 f s samples/s (2) 

Then, a first estimate f 0 is taken as the frequency corresponding to the 
largest FFT output coefficient magnitude. A discriminant which is proportional 
to the frequency error in the initial frequency estimate is computed using 
modified coefficients of the discrete Fourier transform (DFT) with center 
20 frequencies plus one half and minus one half of the FFT bin spacing relative to 
the original frequency estimate f 0 . 

The value of the discriminant is then mapped into the estimate of the 
frequency error in the initial frequency estimate using a mathematically derived 
function. 

25 The estimated frequency error estimate is added to the initial FFT 

estimate to get a next interpolated frequency estimate. 

The process is then repeated, using the next estimate and computing a 
new frequency discriminant to produce a next, more precise, frequency 
estimate. 
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THE FREQUENCY INTERPOLATION DISCRIMINANT 

It is assumed that the signal to noise ratio is sufficiently high such that the 
largest magnitude FFT coefficient corresponds to a frequency closest to the 
5 signal frequency. This assumes that the signal to noise ratio is sufficiently high 
that the probability of the statistical outlier event of a noise only FFT bin 
magnitude being larger that a FFT bin containing both signal and noise is 
neglible. 

10 Define, 

a =Zr[n]exp[-j2^(— -— )n] {7) 



N-l f ] 

P = £ r[n]exp[-j2/r(— + — )n] 

n=0 4 t» *™ 



(6) 



Then the discriminant, or distance metric, of frequency estimation error is 
15 defined as, 

DM =mii«i (9 ) 

where, * = fT f -.^« (10) 
N 

and, 

20 

s N 

k 

For the initial frequency estimate using the FFT, f 0 T s = -ff-and s = 0. 

N 

In the noiseless case, 

25 
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. -1 
-1, s-e = — 
2N 

0, W=0, (11) 

« 1 

1, e-e- — 

2N 



D(e t s) is a monotonically increasing function of e-e . Therefore, each£>(£-,£), 
there is a unique inverse mapping to e - s . Clearly, D(s t e) ) may be used as a 
5 discriminant for fine frequency interpolation between FFT bin centre 
frequencies. 

There exists some functional relationship such that, 

10 ^=V + ^ D(ff '^' (12) 
where, ^( . ) is a monotone increasing function. 
y(.)is called the frequency interpolation function and 
f n is the interpolated frequency estimate. 
The requirement that f< has zero error in the noiseless case is, 

15 iy[D(£j)]=£-£ , for-l^D^l. Therefore, y/' x (£•-<?) = D(s,£) . 

THE FREQUENCY INTERPOLATION FUNCTION 

Assume that, — — <fT <^-+— . 

N 2N 5 N 2N 

20 

Then, — l -~$e-s< — . (13) 
2N 2N 

r(n]=s[n] = Ae j3;rfhT -, n = 0 J l,2,...,N-l (14) 

25 

Without loss of generality, assume that £ = 0. Also assume the noise free 
case. 



The FFT output coefficients are given by, 
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Y(k) = Ai-i^ > k = 0,l,...,N-l 

! ' e (15) 
_ Ac i^ ») sm(^gN) 

sin(^f) 

The discriminant can be expressed as , 

I Y ( k m« +^)|-| Y ( k ma "1 I I S«tl(^ + | " | Sm(^- | 

5 D(*,0) = -2 f- = f f- (16) 

I Y(k raK + ^) I + 1 Y(k mK - i) | | sin(^ + — ) | + 1 sin(^ - — ) | 

After some trigonometric simplification, 
tan(— ) 

D(^,0) = — (17) 
tan(7r£) 

10 This inverse mapping from T>(e, s) to s - e can be obtained as, 

s - £= irfD(£> - tan" 1 [D(* B *)tan(^-)] (1 8) 

Then the frequency estimate, f,,may be obtained, where, 

f T > =^ + £ = iSn^ + ^ (D)= iEm« + 1^^ £) tan (JL)] (19) 

N N N *r 2N /J 

15 

^ = V + ~ tan [ f V^O 1 (20) 

The implication is that a two point (N=2) FFT is sufficient to obtain zero error 
frequency determination in the noiseless case. However, the Cramer-Rao 
20 lower bound is relatively large for N=2 and the SNR threshold is relatively large. 
There is a motivation to use larger N to reduce the rms frequency estimation 
error. However, the implication of larger N is increased computational 
complexity and longer delay time to obtain the transform results. It is desirable 
to obtain low frequency estimation error with the smallest possible N. 

25 




13 

Define, 

r,(D) = ~ ,for-l5DSl. 

For large N or for any N and small |D|, the function y, (P) closely approximates 
5 y(D) . 



f,T, = lim 



N 



',(D) = !im 



N 



|Y(k m 



+|)l + |Y<k. 



-J). 



W 1 



N-»oo 



N-»oo 



(21) 

10 

FREQUENCY ESTIMATION BY ITERATION 

For the case of r[n] consisting of signal plus noise, the noise will cause a 
15 perturbation of D(s,£), and some error in the frequency estimate will result. 

Although an algorithm for the exact frequency determination in the noiseless 
case has been presented, it will be shown that the noise performance of this 
algorithm improves substantially when \s-£\ is close to zero. Since the 
discriminate D(£r,£) can be used to get an interpolated frequency estimate with 

20 less than one half FFT bin size error, it then follows that the algorithm can be 
iterated to move towards zero and |D(s,£)| towards zero. In this way 

the variance of the frequency estimator output can be reduced. 



The iterative algorithm is defined below. 

25 

Define, 

§ r[n]exp[-j2*(^ + -J-)] - |>[n]exp[-j2*(^ - -L)] 

■d(*.o =^ & 5' 2N 
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Define a monotone increasing function of 
' D, such that 0) = 0, *1 ) = and #-1 ) 

and, 

5 . k 

* r y max 

*m = I m A S N 

The iterative algorithm is defined by, 
k 

10 

which implies that s 0 = 0 . 
Then, 

15 f,T.=f 0 T.-^D(f,* 0 ] 
f 1 T,=f 1 T,4tfD(*,ff 1 ] 

and, 

m — > oo 
D^lim (DJ 
m — > oo 

25 The steady state frequency estimate at the end of the iteration is f„ . 

Define, s k =f k T 5 for k=0, 1,2,3 

Then, 





and, the normalized frequency error is, 
-f)T s 



5 

The iteration may be viewed as a convergence to a fixed point of the equation, 

K = g(s^ ) = C - <p{V(* , C )1, formal 

and < £ 0 < 

2N 0 2N 

10 Theorem 1 below is referenced, [21]; 

Let g(x) be a continuous function on [ — — ] and g([a,b])c [ 

2N 2N 

Furthermore, assume that there is a constant 0 < X < 1, with, 
|#(x)-g(y)|<,l|x-y|, forallx,ye[a,b], 
Then, 

x = g(x) has a unique solution x rt in [a,b], 
15 aJso, 

the iteration x n =#(x n .,),for n > 1 will converge to 
x^for any choice of x 0 € [a, J], 
and, 

In the situation under analysis, for fixed { rfn]}* A and s, <p()is a function of e. 



2N'2N 
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and, 

ls(x)-s(y)N|x-y||i.^iM|. 

x-y 

ifo<«»>-*y><2 

x-y 

Then, 

|g(x)-g(y)M x-y|, 

wh ere,A = m ax|l-M^M| forx , ye[ » » j 
x-y 2N 2N 

Using Theorem 1 , it follows that the iteration will always converge to a fixed 
point under the appropriate conditions. 

5 Also, <*(D W ) = 0, and from the properties of {*(.), = 0. 
The fixed point solution, satisfies, 



I £ r[n]exp[-j2;r(^L + J^) n ] | - 1 f>]exp[-j2*(^ - -L)n] | 
D = _Ho f« 2N ft=Q f 5 2N 

10 9 VA f 1 NI f 1 

| £r[n]exp[-j2;r(^ + -i-)n] | + 1 2>]exp[-j2*(f. - — )n] | 

The two previously defined functions ^(D) and ^,(D) fulfil the requirements of 
^(D) and may be used in the iteration to obtain f^. While ^(D) iteration will 
tend to converge more rapidly than ^,(D) iteration, both will yield identical 
15 values of f w . However, evaluation of ^(D)has lower computational complexity 
than evaluation of ^(D). There is performance advantage in using ^(D)when 
the computation is limited to a few iterations. 

PERFORMANCE ANALYSIS 

20 

For the case of signal in additive noise, is a random variable. 
varffJ-J = var[*J= var^D^ )J 
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D(ff, J |/?H«I 

In general, |a| and |P| are both Rician distributed random variables. However, 
tinder high signal to noise ratio, both |a| and |(3| are essentially Gaussian 
5 distributed random variables. 

e rt , which is a random variable, is found by the constraint, 

O.D.-D(r.O-{|{={f{. *rf-f. 
10 ^will be perturbed by the noise component in D,. Even though D^is 
constrained to be zero, the constraint and noise induce randomness in e^. The 
noise perturbation in 
D M induces the perturbation in . 

15 The approach taken is the computation of the variance of Dfrom the point of 
view of the creation of D from noisy observations and then to find the 
corresponding perturbation of - s . 

For high signal to noise ratios, Dfo^Jmay be represented by the first three 
20 terms of the Taylor Series expansion about, ix a and p fi , which are the means 

of |a| and | /? |, respectively. 

\(3\ + \a\ Mfi+M a (M fi +M a ) (M 0 +fa) 

25 Assuming £„, = £ 
A 

Ma =M P = 



sinf— ) 
2N 



and for high signal to noise ratio, 

-1 
2 



30 Then, 
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E[D]=0, 

varfD] = + 

It then follows that, 

5 

Nsin 2 (^) 

var[D] = 2fiL, 

J 4(SNR) 

where, SNR = ^7- 
cr 

10 = — tan* } [D tan(— )] 

n 2N 

Then, for high signal to noise ratio, the normalized frequency error may be 
computed. The largest part of the probability density function of D is in the 
region of where the atan(x) * x . Therefore, 

15 

E[atan 2 [Dtan(^)] <r£tan 2 (-^) Nsin 2 Atan 2 (^~) 
<7;-var(*-*)- - 2 ? - 4 (SNR)^ 



PERFORMANCE COMPARISON TO THE CRAMER-RAO LOWER BOUND 

20 

The Cramer-Rao Lower Bound on the variance of the frequency error of any 
unbiased frequency estimator is given by, 



(2*r) 2 N(N 2 -l)(SNR) 
25 

The performance of the DFT based estimator may be compared to the Cramer- 
Rao Lower Bound. 
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a] ^(N^)sin^)tan^) 

^CRLB 6 

For high SNR and large N, the performance frequency estimation variance is 

10log 10 ( — ) = 0.063282577 dB .above the Cramer-Rao Lower Bound. 
96 

2 

5 Figure 5 shows 2 * in dB verses N t where N is the length of the FFT. 

&CRLB 

JUSTIFICATION FOR ITERATION TOWARD A ZERO VALUE OF THE 

DISCRIMINANT 

10 

The reason for the performance improvement of the proposed class of 
algorithms relative to prior algorithms is the first frequency interpolation allows 
the computation of two DFT coefficients, which are !4 DFT bin spacing above 
the first interpolated frequency and V* DFT bin space below the first 

15 interpolated frequency. While the first interpolation may still have significant 
error, which is dependent on the relationship of the true frequency relative to 
the FFT coefficient frequencies, the error discriminant evaluated for the first 
interpolated frequency will have a value close to zero. The variance of the 
frequency error is relatively low in the region of small values of the frequency 

20 discriminant. Therefore, the second interpolated frequency will have small error 
variance. There is significant noise performance advantage in using the first 
interpolation to allow a low error variance second interpolation. The 
interpolation may be iterated, with diminishing improvements of estimation 
accuracy, until convergence to a fixed point solution is obtained. Figure 6 

25 shows the variance of the normalised estimator frequency error estimate vs the 
frequency error for the first interpolation. For the figure, N=64 and the signal to 
noise ratio is 6 dB. There is a very sharp reduction the rms error of the 
frequency estimator in the region of the frequency being close to the center 
frequency of the frequency discriminator. This indicates that tremendous 

30 improvement in performance obtained by iteration. 



ALGORITHM SIMPLIFICIATONS RESULTING IN LARGE REDUCTIONS IN 
COMPUTATIONAL COMPLEXITY 
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Simulation results have verified a single FFT and two iterations involving the 
computation of the discriminant function D(£-,^ 0 )andD(£-,^)are sufficient to 

obtain rms frequency error performance close to the computation of f m . The 

5 first iteration moves the discriminant towards zero and decreases of | e-e x |. 

The estimate resulting from the second iteration therefore results in small error 
variance of € - s 2 . 

The algorithm has complexity of 0(N log 2 (N) +0(8N)=0[N log 2 (N)] 

10 

The algorithm has the same order of complexity as the original FFT for 
performance, which is very close to the Cramer-Rao Lower Bound. 

The algorithm has the same order of complexity as the original FFT for 
15 performance, which is very close to the Cramer-Rao Lower Bound. 

Iteration with ^,(D) =D yield results very close to ^(Z))and saves considerable 
computational complexity. The fixed point solution for the two cases is 
identical. 

20 

Since the algorithm comes very close to convergence with two iterations, two 
iterations are sufficient for most applications. The performance improvement 
with additional iterations is small. 

25 

OTHER FREQUENCY DISCRIMINANTS WITH THE SAME NOISE 
PERFORMANCE 

There are a number of discriminants, which have the same performance, when 
30 used iteratively to obtain the fixed point solution, as the previously introduced 
discriminants. The noise performance is identical, for iteration, because the 
fixed point solution is identical. 

This class of discriminants includes functional forms, 

35 





5 



and in particular, D = 



and, 



D = Re[- 



/3 + a 



T ] , where * denotes complex conjugate. 



for/>0. , 



WPf-\a\ 2 
2\p?+\a? 
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And where Re[.] is the real part. 



10 



SIMULATION 



Figures 7 to 12 show the rms frequency error performance of the algorithm vs 
SNR in dB, for N=2,4, 16,64,246, and 1024, respectively. Both the cases of one 
15 interpolation and two iterative interpolations are shown. The two interpolation 
case is essentially achieves the performance of the infinite interpolation case. 



In summary, a new, low computational complexity, class of algorithms, which 
interpolates the result of a FFT, has been presented for the precise estimation 

20 of frequency of a complex exponential function in additive Gaussian noise. The 
performance of the algorithm, above the threshold in additive Gaussian noise, 
is about 0.06 dB above the Cramer-Rao lower bound. The algorithm is ideally 
suited to be utilized in a number of communications, signal processing and 
biomedical applications. The algorithm also has ideal characteristics for digital 

25 signal processor implementation. 



Industrial Application 

There are a large number of applications for this invention, including: 
30 - Rapid frequency initialisation of a phase lock loop for rapid signal acquisition; 
- Radar processing for precision radial velocity and radial acceleration target 
measurements; 
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- Sonar processing for precision radial velocity and radial acceleration target 
measurements; 

- Satellite orbit determination; 

- Missile trajectory determination; 

5 - Ultra sound imaging Doppler measurements for blood and other biological 
fluid velocity measurements; 

- Ultra sound imaging for tomography processing involving Doppler shift 
measurements; 

- Coherent carrier tracking for coherent digital demodulators - large frequency 
10 acquisition range and rapid signal acquisition; 

- Noncoherent digital communication system frequency tracking - large 
frequency acquisition range and rapid signal acquisition; 

- Frequency estimation for electronic test equipment displays including 
frequency meters, oscilloscopes, spectrum analyzers and network analyzers; 

15 - Ultra low distortion, ultra high performance FM demodulator, and 

- Generalised software modules in Matlab and other commercial software 
packages. 

It will be appreciated by persons skilled in the art that numerous 
20 variations and/or modifications may be made to the invention as shown in the 
specific embodiments without departing from the spirit or scope of the invention 
as broadly described. The present embodiments are, therefore, to be 
considered in all respects as illustrative and not restrictive. 

25 Dated this nineteenth day of July 2002 

University of Technology, Sydney 
Patent Attorneys for the Applicant: 



30 



F B RICE & CO 
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